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The Glauber method is extensively used to describe the motion of a hadronic 
projectile in interaction with the surrounding nuclear medium. One of the main 
approximations consists in the linearization of the wave equation for the interacting 
particle. We have studied the consequences of such an assumption in the case of the 
12 C(e, e'p) 11 B* reaction at high proton momenta by comparing the results with the 
predictions obtained when all the ingredients of the calculation are unchanged but 
the second-order differential equation for the scattered wave, which is solved exactly 
for each partial wave up to a maximum of 120 spherical harmonics. We find that 
the Glauber cross section is always larger by a factor 10 -j- 20%, even at vanishing 
missing momenta. We give a quantum-mechanical explanation of this discrepancy. 
Nevertheless, a good correlation is found between the two predictions as functions 
of the missing momentum, especially in parallel kinematics. 
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finite nuclei is receiving more and more interest in relation with the experiments planned 
at CEBAF, where the proton momentum can be larger than 1 2 GeV/c. At large 
missing momenta of the recoil the details of the short-range nucleon-nucleon interaction 
are expected to show up 0, but the few experimental data available || still prevent from 
putting stringent constraints on the theoretical models. On the other side, at moderate 
missing momenta an accuracy within 10% is required to unambigously identify exotic 
effects like Colour Transparency, if any [|J. 

The recent NE18 experiment |J has shown that calculations based on the standard 
Glauber approximation overestimate the cross section even at vanishing missing momenta, 
unless some phenomenological corrections are introduced. Various suggestions have been 
made both on the way of analyzing the NE18 data || [7j and on how to improve the 
Glauber method ||. Starting from a different point of view, we show that the latter 
is affected by an intrinsic systematical error, which can be relevant at the kinematics 



1 



explored at CEBAF but seems anyway to be predictable and perhaps can be corrected 
for. 

Traditionally, the Glauber method has been extensively used in the analysis of the 



proton-nucleus scattering data || [UJ, where it gives good, or at least some of the best 
available, fits to the scattering distributions. However, the generalization to the (e,e'p) 
reaction is not straightforward. In fact, even if the hadronic content of the final state is 
similar, the inelastic nature of the electromagnetic proton emission makes the kinematical 
situation rather different. 

The basic ingredient to describe the knockout process for an exclusive scattering in 



the framework of DWIA, is the scattering amplitude |Tl], [12 



J£(ff) = Jdrdae^'^^^J^q^a^^a) , (1) 

where q is the momentum transfer, x is the distorted wave function of the knockout 
nucleon and is the single-particle bound state wave function depending on the quantum 
numbers a of the hole. For sake of simplicity, in this work we have focussed on the 
longitudinal component of the current operator, retaining just the leading order 0(1) in 
the nonrelativistic expansion and neglecting the nucleon form factor, i.e. substituting J M 
with the identity operator. The cross section becomes, therefore, proportional to 

Jdrdae^-^^a^^atf^S^q) , (2) 



which can be identified as the "distorted" spectral density [fnj at the missing energy 
corresponding to the knockout hole a. 

The distorted wave function x is solution of the Schrodinger equation 

("2^ V2 + V ) X = EcmX ' (3) 

where m is the reduced mass of the nucleon in interaction with the residual nucleus, E cm 
is its kinetic energy in the cm system and V in principle includes a nonlocal energy- 
dependent optical potential effectively describing the residual interaction. The incoming- 
wave boundary conditions to be imposed on the solution of eq. (D) correspond to the 
well known feature of the quantum-mechanical problem for a wave crossing a potential 
barrier, where the asymptotic stationary conditions require the incoming unitary flux to 
be splitted in a reflected one and in an outgoing one. 

The proper way of solving eq. (|3|) is to expand x^ in partial waves, solve the second- 
order differential equation wave by wave, sum the solutions up to a certain maximum 
angular momentum L max which satisfies an ad hoc convergency criterion for eq. ([I]). 
From now on, this procedure will be referred to as method I. So far, it has been applied 
to proton momenta below 0.3 GeV/c [E|, where convergency is reached for L max < 50. 



In this work we have extended it to larger proton momenta, i.e. p ~ 1 4- 2 GeV/c. A 
maximum L max = 120 has been used and we have checked that the results become already 



stable with L max = 100. 
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At high energies the Glauber method JIJ offers an alternative way (from now on method 
II) of solving eq. (|3|) by linearizing it along the propagation axis z: 
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where b is the impact parameter describing the degrees of freedom transverse to the motion 
of the struck particle with momentum p. With this approximation eq. (|3[) becomes 

The standard boundary condition applied requires that asymptotically ^ — > 1, which 
corresponds to an incoming unitary flux of plane waves. In contrast to method I no 
reflected flux is taken into account that is created by the potential of the nucleon-nucleon 
interaction. Because of the different phases, there is a destructive interference between the 
reflected and the incoming waves, which eq. (|7|) neglects thus leading to an overestimation 
of eq. In particular, the approximation of eq. @, upon which the eq. (0) is based, 
can be expected to produce most of the discrepancy at very small angles 7 between q 
and p, typically in parallel kinematics. On the other side, at large 7 and large transverse 
missing momenta the approximation of eq. ([5]) is expected to play the major role. 

A direct comparison with the exact solutions of eq. (|3|), summed up to L max , is 
needed to get a quantitative answer out of these qualitative considerations. The need 
for considering corrections beyond the eikonal approximation of the Glauber method is a 



well known problem in high-energy elastic proton scattering on nuclei |]10, [T4|, 15 . One 



basic outcome of these studies is that competing effects tend to cancel each other at small 



deflection angles pM , while can give important corrections at larger angles |[L5| . This fact 
suggests that also in (e,e'p) reactions it could be interesting to compare deviations from 
the eikonality with other kinds of corrections to the Glauber model. This is the aim of a 
more general work, whose preliminary results we present in this communication. 

Here, we have considered the 12 C(e,e'p) 11 B* reaction both in parallel and perpendic- 
ular kinematics for the proton momentum p up to 2 GeV/c and the momentum transfer 
q = 1.4 GeV/c. When not explicitely mentioned, the bound state ty Q must be understood 
as the solution of the Woods-Saxon potential of Comfort and Karp |16| with the quantum 



numbers a of the s wave. As an exploratory calculation, we have neglected the contri- 
bution to V in eq. (0) coming from the Coulomb potential to avoid numerical problems 
related to the high angular momenta required. Therefore, in proper terms the results 
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presented here refer to the (e,e'n) reaction. V is an optical potential of the form 



V(r) = (U + iW)- — 

1 + e— 

= (U + iW)p(r) , (8) 

with R = 1.2 x A 1 / 3 fm and a = 0.5 fm. The nuclear density p(r) defined in eq. (§) is 
normalized such that p(0) = 1. 

No phenomenological optical potential is available at the energies here considered. 
Therefore, the parameters U, W can only be guessed. Two guidelines can help in this 
case. On one side, according to the Glauber model the imaginary part should scale as 
W ~ p/10 MeV, while U/W should equal the ratio between the real and the imaginary 
parts of the average proton-nucleon forward-scattering amplitude, which is expected to be 



< 0.5 in the considered kinematics |l^j. On the other side, assuming that the absorption 
due to the final state interactions, observed in the NE18 experiment in the context of a 
semi-inclusive (e,e'p), is reasonable also for a completely exclusive knockout, then U, W 
should be better taken about half of the previous values. We have considered several 
choices, including the unphysical but interesting cases of a completely real (W = 0) or 
imaginary (U = 0) potential. 

By comparing the predictions of methods I and II for eq. (0) at p = q we find that in 
the range 0.6 < p < 2 GeV/c the results are p-independent if U, W are linear functions 
of p. This is a well known feature of the Glauber method at high energy proton-nucleus 
scattering. Hence, the previous findings suggest that, at least for small missing momenta, 
starting from p ~ 0.6 GeV/c on we are already in a "Glauber regime". Therefore, in the 
following whenever the kinematics is parallel, we employ an optical potential of the form 
(U + iW)p{r)p/p , where p a = 1.4 GeV/c. 

Figs. 1 and 2 show the distorted spectral density of eq. (0) in parallel kinematics for 
missing momenta in the range < p m = p — q < 600 MeV/c. The dashed line represents 
the result with no final state interaction (PWIA), which is of course identical in both 
methods. The solid and the dotted lines are the results of method I and II, respectively. 
The two curves are rather well correlated in all the range of p m here explored. But a 
systematic discrepancy is evident, even at p m ~ 0, which is roughly proportional only to 
the imaginary part of the optical potential, as it can be realized by inspecting curves in 
fig. 1 labeled by (a) (U = 100 MeV, W = 0) and (b) (U = 0, W = 100 MeV), or in fig. 2 
by (a) (U = 20, W = 50 MeV) and (b) (U = 50, W = 150 MeV). This fact suggests that 
the origin of the discrepancy should be ascribed to the linearization of eq. (^) implied by 
the approximation of eq. (||). 

It is interesting also to explore the role of the Final State Interactions (FSI), which 
seem to have a nontrivial structure especially at high missing momenta in perpendicular 
kinematics. To accomplish this, we rewrite eq. (0) as 

S%(q) ~ |PWIA + FSI| 2 

= |PWIA| 2 + | FS 1 1 2 + 2Re(PWIA • FSI*) (9) 
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and we compute the single contributions separately in the framework of method I. In fig. 
3 the |PWIA| 2 , | FSI | 2 and the total contribution are shown as functions of the angle 7 
between q and p using a strong unrealistic optical potential (U = 50, W = 150 MeV) to 
emphasize the role of FSI. At small angles the PWIA is the dominant contribution. The 
interference is negative and is responsible for the familiar damping observed in all electron 
scattering data. Near 7 ~ 15° the angular distribution has a dip (partially filled if the real 
part U of the optical potential is non vanishing), because the interference is roughly equal 
and opposite to the sum of |PWIA| 2 and |FSI| 2 . At large angles the distribution basically 
coincides with |FSI| 2 . 

To better specify this last issue let's consider the fig. 4, where the angular distribution 
is plotted in the same conditions as in fig. 3 but substituting inside \1> Q the shell- 
model wave function of the Woods-Saxon type with a harmonic oscillator with the same 
quantum numbers. In this case the dashed line, i.e. the PWIA contribution, falls down 
monotonously and quickly just after the Fermi momentum (221 MeV/c at 7 ~ 10°). Any 
diffractive pattern must therefore be ascribed to the FSI term. This remark, less obvious 
in the case of fig. 3 because of the richer structure of PWIA, is a common feature of 
both approaches, the method I (solid line) and the method II (dotted line). They appear 
correlated similarly to the case of the parallel kinematics. 

The oscillatory diffractive trend of at large angles is reminiscent of the angular 
distribution for proton- nucleus scattering |J. Thus, the natural interpretation is that 
the ejected proton is testing coherently the residual nucleus. This is possible only in 
a completely exclusive reaction, i.e. where the residual nucleus doesn't fragment. Any 
measurement of an energy-integrated distribution, i.e. of a semi-inclusive (e,e'p) reac- 
tion H [7], |18|, [19j , looses this information on the structure of the A — 1 system and tests 



just the average behaviour of the single ejected nucleon. 

We have critically analyzed some aspects of the Glauber approach to the description of 
the propagation of a particle through the nuclear medium in the case of the 12 C(e,e'p) n B* 
reaction. The basic approximation of linearizing the wave equation for the projectile has 
been shown to be responsible for a systematic discrepancy with respect to the results 
obtained from its proper solution. This fact could be relevant at the CEBAF kinematics, 
particularly at small missing momenta. At very large values of the transverse missing 
momenta the FSI give the dominant contribution to the cross section and show a diffractive 
pattern, which can naturally be interpreted as a coherent diffractive scattering between the 
ejected proton and the residual nucleus. To verify this prediction a completely exclusive 
(e,e'p) reaction is needed, where the residual nucleus is in a well specified state. 

We would like to thank prof. S. Boffi for many discussions and for his continuous 



interest in this work, and dr. F. Cannata |20| for having originally suggested the idea 
upon which this work is based. 
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Captions 



Fig. 1 - The distorted spectral density as a function of the proton momentum p for 
the 12 C(e, e'p) n B*iy 2 reaction in parallel kinematics at the momentum transfer q = 7fm _1 . 
The dashed line is the result of the PWIA. The solid and the dotted lines are obtained 
in the framework of method I and II, respectively (see text). The curves labelled by (a) 
are produced with U = 100, W = MeV/c, where U, W are the real and the imaginary 
parts of the optical potential, respectively. The curves labelled by (b) are produced with 
U = 0, W = 100 MeV/c. 

Fig. 2 - The same as in fig. 1, but the curves labelled by (a) are produced with U = 
20, W = 50 MeV/c, while the curves labelled by (b) are produced with U = 50, W = 150 
MeV/c. 

Fig. 3 - The as a function of the angle 7 between q and p for the 12 C(e, e'p) 11 B* 1/ , 2 
reaction in perpendicular kinematics at q — 7fm -1 . The dashed line is the PWIA result, 
the dotted line is the pure contribution of the FSI (when the PWIA contribution is 
subtracted), the solid line is the coherent sum of the two. All the curves are obtained in 
the framework of method I with U — 50, W — 150 MeV/c. 

Fig. 4 - The same as in fig. 3 but with the bound state described by a harmonic 
oscillator. The dashed line is the PWIA result, the solid and dotted lines are the total 
result for method I and II, respectively. 
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